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Filtering  by  means  of  an  analog  computer  to  separate  an  unwanted  random  signal  from  a 
lower  frequency  signal  is  discussed.  The  spectral  concept  of  the  filtering  problem  is  men- 
tioned and  an  example  of  an  analog  computer  filter  circuit  is  given. 


INTRODUCTION 


The  smoothing  of  electrical  signal  data  by  filtering  on  an  analog  computer  is  not  a particularly  new 
procedure.  There  are,  of  course,  more  sophisticated  analog  methods  of  extracting  a desired  signal  which 
is  intermingled  with  other  signals.  However,  we  shall  briefly  discuss  the  subject  from  the  standpoint  of 
filtering  to  separate  an  unwanted  random  signal  from  a more  slowly  fluctuating  signal.  We  hope  that  this 
paper  will  have  added  benefit  for  the  person  who  has  not  previously  encountered  references  on  the  general 
technique. 

The  filtering  method  presented  was  used  by  the  second  author  for  smoothing  the  output  signal  from  a 
force  transducer  in  coefficient  of  friction  measurements  for  agricultural  product  materials  and  by  the 
first  author  to  treat  densitometer  signals  in  the  determination  of  background  densities  on  spectrographic 
films. 

The  reader  will  find  the  references  by  Korn  and  Korn  ^ and  Jackson  ® helpful  in  regard  to  the  simulation 
procedures  used. 


GENERAL  THEORY 


In  the  case  under  discussion  it  is  desired  that  the  output  signal  should  have  a narrow,  low-frequency 
spectral  density  function  T^(cu).  For  an  input  signal  spectral  density  <I>(a>),  the  hnear  filter  of  character- 
istic 0(ct>)  will  operate  on  0(o))  in  the  manner 

^(o>)  = 0(&>)<I>(a>).  (1) 

We  now  consider  the  filter  analogous  to  the  system 

V(t)+aV(t)  = U(t).  (2) 


^ Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles,  Agricultural  Engineering  Research  Divisions,  ARS,  USDA,  and 
the  Department  of  Agricultural  Engineering,  Ohio  Agricultural  Experiment  Station,  Wooster,  Ohio,  respectively. 

^ G.  A.  Korn  and  T.  M.  Korn.  Electronic  Analog  Computers.  2d.  ed.  New  York.  1956. 

^ A.  S.  Jackson.  Analog  Computation.  New  York.  1960. 


1 


(3) 


In  the  manner  of  Batchelor,^  we  write 

V{t)=  f?^‘^dZ(co);U(t) 

—00 

and  define 


hm 

d(i) 


dZ{(x)i)dZ*{(x)2) , 

d(ii 


0; 


lim 

doi- 


dWMdW*{0i2)_r. 

0 

(0)l  7=  (O2) 


= dZ(w)dZ*(co)  0(0,)  = ^"'  dmo>)dlP'*M 

d<o-*0  ^ ’ dw-^0 

Equations  (3),  (4a),  and  (4b),  when  apphed  to  (2),  produce  the  result 


(ift>  + A.)dZ(a>)  = dW{(xi), 


and  further. 


Hence,  we  finally  obtain 


(o)2  + \^)dZ{oi)dZ^{oi)  = dW{(x))dW*{(xi). 


Comparison  of  (1)  and  (7)  shows  that 


Oioj) 


1 


0>2  + X2 

represents  a low-pass  filter  which  will  give  a ^'^(aj)  of  the  desired  type. 


THE  FILTER  CIRCUIT 


(4a) 

(4b) 

(5) 

(6) 

(7) 

(8) 


In  analogy,  suppose  we  wish  to  examine  the  velocity  of  a sphere  of  mass  m subjected  to  an  arbitrary 
force  u{t)  and  immersed  in  a viscous  fluid.  We  presume  that  Stokes’  law  dissipation  exists,  i.e.,  the  vis- 
cous force  Fv  will  be  given  by 

Fv — — K{dxl  dt).  (9) 

If  random  fluctuations  £U‘e  superimposed  on  u{t),  it  will  be  obviously  difficult  to  measure  the  velocity. 
However,  if  u{t)  remains  unchanged,  and  we  place  the  sphere  in  a more  viscous  fluid,  the  random  fluctua- 
tions in  the  velocity  will  be  more  heavily  damped  in  accordance  with  the  differential  equation. 


(Fx  Y.dx  1 . . 

r — — 

dt^  m dt  m 


For  the  velocity  v = dxjdt,  (10)  becomes 


dv  K _1  .. 

-r-\ — v= — u(t). 
dt  m m 


(10) 


(11) 


■*  G.  K.  Batchelor.  The  Theory  of  Homogeneous  Turbulence.  London.  1953. 
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We  note  that  (11)  is  of  the  same  form  as  (2).  Therefore,  placement  of  the  sphere  in  the  more  viscous 
fluid  would  he  analogous  to  the  use  of  a low-pass  filter.  Taking  the  Laplace  transform  of  (11)  and  assuming 
initial  conditions  are  zero,  we  find  that 


y(s) 


K(f+.) 


U(s). 


(12) 


The  first  factor  on  the  right-hand  side  of  (12)  is  the  type  of  transfer  function  needed  to  simulate  this 
low-pass  filter.  Figure  1 shows  the  appropriate  operational  amplifier  connection,  for  which  the  output 
^0  in  terms  of  the  input  Xi  is  ® 


Rq Y 

«(«oCs+l) 


(13) 


With  Ro  constant,  the  capacitance  C is  directly  proportional  to  the  time  constant  of  the  circuit,  which  in 
turn  governs  the  frequency-cutoff  point  of  the  filter. 

The  circuit  in  figure  1 does  not  have  the  sharp  cutoff  characteristics  of  some  more  complex  filters. 
However,  it  was  very  adequate  for  this  work  and  there  was  good  recorder-driving  capabihty. 


APPENDIX 


An  example  of  the  use  of  the  filter  circuit  in  Figure  1 was  the  separation  of  unwanted  higher  frequen- 
cies from  the  desired  low  frequency  signal  of  a friction  coefficient  force  transducer.  The  circuit  was 
implemented  on  a Boeing  Model  6631  analog  computer.® 


1^0 


® See  p.  415  of  footnote  2. 

®The  mention  of  trade  products  does  not  imply  that  they  are  recommended  or  endorsed  by  the  Department  of  Agriculture  over 
similar  products  of  other  companies  not  mentioned.  Trade  names  are  used  here  for  convenience  in  reference  only. 
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The  value  of  the  resistors,  R and  Ro,  were  set  at  1 megohm  each  since  these  values  were  readily 
available  within  the  operational  amplifier.  The  value  of  the  capacitor,  C,  was  adjusted  from  0.1  to  0.9 
mfd  to  determine  the  optimum  degree  of  filtering.  The  value  of  0.5  mfd  gave  the  best  balance  between 
the  response  and  smoothing.  As  shown  in  figure  2,  the  values  of  0.9  and  0.7  mfd  for  C were  excessive 
since  the  filtered  trace  did  not  reach  its  maximum  before  the  end  of  the  test  run,  whereas  the  value  of  0.5 
mfd  did  reach  a maximum  soon  after  the  start  of  the  test  run. 


0.5  afd. 


Figure  2. — Filtered  and  unfiltered  transducer  output  as  obtained  with  various  values  of  the  capacitor,  C. 

filtered  output  is  on  the  right  and  the  unfiltered  output  is  on  the  left. 


In  each  set  of  curves  the 
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